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^ ' Abstract: We prove that the time derivative of the solution for the obstacle 

problem related to the Evolutionary p-Laplace Equation exists in Sobolev's sense, 
provided that the given obstacle is smooth enough. We keep p > 2. 
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1 Introduction 



\^ • The celebrated Evolutionary p-Laplace Equation is much studied and the reg- 

GN . ularity theory for the solutions is almost complete. We refer to the book [dB] 

pq . about this fascinating equation. In general, the corresponding subsolutions 

(X '■ and supersolutions do not possess that much regularity, they are semicon- 

O I tinuous. We are interested in a special kind of weak supersolutions of the 



Evolutionary p-Laplace Equation, namely the solutions of an obstacle prob- 
lem. In the presence of a smooth obstacle the regularity improves a lot. 
Given a function ip = il){x,t) in a bounded domain VLt = fi x (0,T), where 
VL C R*^, we consider all functions v such that 

ov 

— >V-{\VvY'-'^Vv) and v > i) in VLt- 

The function ip acts as an obstacle. The smallest admissible v is the solution 
of the obstacle problem. (This makes sense because a comparison principle is 
valid.) However, the above description was only formal. We will instead use 
Definition 1 below, which is more adequate since it comes with a variational 
inequality. — We will restrict ourselves to the case p > 2, the so-called slow 
diffusion case. 




It is an established fact that if the obstacle ip is smooth enough, the 
solution to the obstacle problem inherits some regularity. Our objective is 
the time derivative Ut of the solution u, which a priori is only known to be 
a distribution. Our main result Theorem 2 states that, if ip has continuous 
second derivatives, then the time derivative Ut exists in Sobolev's sense and 
it belongs to the space L^i^c (^t)- A formula is given for the derivative. 
The most laborious part of the proof is to show that Ap-u = V ■ (| Vf |*'~^Vf ) 
is a function so that the rule 

T 

ip\u dxdt = - I I { I Vm|^" Vm, Vv?) dx dt 
Q Jo Jn 

with test functions applies. The equation has first to be regularized, keeping 
the obstacle unaffected, and then difference quotients are used. The test 
functions in [LI] can be adjusted to work here. 

An important feature, typical for obstacle problems, is that in the open 
set T = {u > ip} where the obstacle does not hinder, u is, actually, a solution 
to the differential equation. Thus in T the equation Ut = ApU holds in the 
weak sense. The boundary of the coincidence set E = {u = ip} is crucial. 
This enables us to get an identity for the integral JJ mpt dx dt, from which 
one can deduce the existence of the time derivative sought for. The special 
case with no obstacle present was treated in [L2]. — See also [BDM] for some 
general comments valid for "irregular" obstacles. 

To this we may add a curious fact valid for ip G C'^{Q.t)- At all points in the 
coincidence set H the obstacle satisfies the inequality 

dip » , 

Thus a point at which -M- < lS.pip cannot belong to the coincidence set. This piece 
of information follows from the characterization of continuous supersolutions as 
viscosity supersolutions, cf. [JLM]. Then ip itself can do as a test function for the 
pointwise testing required in the theory of viscosity solutions. (The reader may 
consult [K] for some basic concepts.) — We will not need this observation. 

It is likely that the time derivative belongs to the space L'^^^iVLT), but an 
eventual proof of this improvement would require much stronger regularity 
considerations for Vm. We have kept p > 2, but one can expect a counterpart 
to Theorem 2 valid in the extended range p > 2n/{n + 2). The difficulty 
about further generalizations with ApU replaced by some operator div Ap is 



the following. It is absolutely necessary that the solutions of the differential 
equation 

— = dwAp{x,t,u,\/u) 

enjoy the property of having a time derivative themselves, in order that the 
corresponding results could be extended to the related obstacle problem. 
This considerably restricts the possibilities. 



2 Preliminaries 

Let Qhe a bounded domain in the n-dimensional space R" having a Lipschitz 
regular boundary. Suppose that a function ip = iIj{x, t) is given in the closure 
of the space-time cylinder Qt = f2 x (0, T). The function ip acts as an obstacle 
so that the admissible functions are forced to lie above ip in Qt- We make 

the 

Assumption: ip e C{Ih) n W^^P{nT). 

For simplicity the obstacle ip also determines the values of the admissible 
functions on the parabolic boundary 

Tt = ^ X {0} U dn X [0,T]. 

The class of admissible functions is 

J^^ = {v e LP{0,T;W^'''{n))\v e C{Il^), v>-^mnT, v = ip on Tt}. 

— We keep p > 2. 

Definition 1 We say that the function u G J-"^ is the solution to the obstacle 
problem, if the inequality 

I I (dV^I^^-^y^, V(0 - U)) + (0 - m) ^yxdt 

>i fi(Pix,T)-uix,T)fdx (1) 

•/Si 

holds for all smooth functions G J-"^. 



The solution exists and is unique, cf.[AL] and [C]. See also [KKS]. It is also 
a supersolution of the equation Ut > ApU, i.e., 



r f({\vu 

Jo Jn ^ 



P-^Vu,Vip)-u^)dxdt>0 (2) 



for all non- negative ip G C^{Qt)- 

Notice that nothing is assumed about the time derivative Uf. Our main 
result is the theorem below. 

Theorem 2 The time derivative Ut of the solution u to the obstacle problem 
exists in the Sobolev sense and Ut G -^fo^ (^t)- It is the function 

^t in S 

ApU in Qt \ 2 

where S = {u = ip} denotes the coincidence set. 

In order to avoid the difficulty with the "forbidden" time derivative Ut in 
the proof, we have to regularize the equation, keeping the obstacle unchanged. 
We replace |VmP~^Vm by 

\^ 



to obtain an equation which does not degenerate as Vm = 0. 
Lemma 3 There is a unique u'^ G J-"^ such that 



'0 Jn 




(^{\Vu'\^ + £2)'^VM^ V(0 - u')) + (0 - u') ^yxdt 

>l f {<P{x,T)-u'{x,T)fdx (3) 

^ Jn 

for all smooth functions in the class T^. In the open set {-u^ > ■0} the 
function u^ is a solution of the equation 



In the case e ^ we have u" G C°°{^t) and ^ G L'^{Q:t)- 

4 



Proof: The existence can be extracted from the proof of [AL, Theorem 
3.2]. The regularity for the nondegenerate case e 7^ is according to the 
standard parabohc theory described in the celebrated book [LSU]. The proof 
of the Holder continuity for the degenerate case e = is in [C] . 

When £ 7^ 0, we can rewrite equation (3) in the more convenient form 




( (I Vm"|^ + e^)~Vu', Vr])+r]^-)dxdt>0 (4) 

V ^ ^ y Ot / 

Ae(x,i) 

valid for all test functions 7] such that 1] > ip — -u"^ in Qt and rj = on Ft- 
We may even use any continuous r] G L^i^O, T; Wq'^^Q)) with ri{x, 0) = 0. 

In order to proceed to the limit under the integral sign in the forthcoming 
equations we need the convergence result below, where u denotes the solution 
to the original obstacle problem, the one with e = 0. 

Lemma 4 

lim f I (\u' - m|p + \Vu' - Vu\A dx dt = 0. (5) 

Proof: It was established in [KL, Lemma 3.2] that 

lim / / \Vu'-Vu\Pdxdt = 0, (6) 



'^^o JO Jn 



but the strong convergence of the functions themselves requires, as it were, an 
extra compactness argument. Since m^ is a weak supersolution, there exists 
a Radon measure /i^ such that 

/ I U\Vu'\'^ + e^)^Vu',Vip)-u'-^)dxdt= I ipdfXe 

for all functions ip G C^{Qt), whether positive or not. This is a consequence 
of Riesz's Representation Theorem, cf. [EG, 1.8]. See [KLP] for details. 

Given a regular open set (for example a polyhedron) U GG Qt, we have 
to verify that 



with a bound independent oi e, < e < 1. Then the lemma follows as in 
[KLP, pp. 720-721]. (There [S] was used.) To this end, choose a test function 
V? G C^{Qt) such that < ip < 1 and ip = 1 inU. A rough estimation yields 







,e 11 

OO' 



By the maximum principle Hm^Hoo < ll'^lloo and ||V-u'^||^p/^ n is uniformly 
bounded, since the gradients converge strongly. This yields the bound Mu. 

D 



3 The gradient estimate 

In order to prove that ApU is a function, u denoting the solution to the 
obstacle problem, we show that the function F = |V'u|*-^~^'*'^V'U, where the 
usual power p — 2 has been replaced by (p — 2)/2, is in a suitable first order 
Sobolev x-space. This will immediately imply the desired result. At a first 
reading one had better to assume that the obstacle ip is as smooth as one 
pleases, say of class C'^{VLt)- Actually, only the Sobolev derivatives ipxix- and 
ipx^ are needed, while iptt does not appear at all. We recall our assumption 
■il) G CiyiT) n W'^'^iyiT) and use the abbreviation 

Under these assumptions about the obstacle i/j = iIj{x, t) we have the following 
result. 

Theorem 5 For the solution u to the obstacle problem, the derivative DF 
of 

F= \Vu\'^~Vu 




exists in Soholev's sense and belongs to L^f,^^ (^t)- The estimate 

' C\DF\^dxdt<C f f {C+\^Cn\^<dxdt 
'0 JQ Jo Jn 

+C f f C\^u\^dxdt + C f /" |VCr|VV^|Prfxrft 
Jo Jn Jo Jn 

+C [ [ Ci\D^^f + \^^t\^) dx dt + C [ C\^^{x,T)\^dx 
Jo Jn Jn 

holds for each non-negative test function ( = C{x) in C^{Q); and C = C{p). 

Proof: The proof is based on the regularized obstacle problem and equa- 
tion (4), where we abbreviate 

A,{x,t) = {\Vu'\^ + e^)^Vu'. 

We denote its solution by u, suppressing the index e. Thus u means -u^, to 
begin with. Given (, the variable x is given a small increment h so that the 
test function 

Tj = ip{x, t) — u{x, t) + C,{xY[u{x + h,t) — ip{x + h, t)] 



= ({xY [u{x + h,t)- u{x, t)] - ({xf [ipix + h,t)- ■^{x, t)] 

-(i-((xr)Hx,t)-^{x,t)] 

is admissible in the regularized equation 



j j (^{A,{x,t),Vv)+V^)dxdt > 0. (7) 



Inserting the test function, we obtain 



J J ((A,(x,t), V(C^Ahix)) + CAhM^) dxdt 

-J J ((A,(x, t), V(C^Ah^)) + C'AhV^^) dx dt 

> f f({A,{x,t),V{{l-C{xrMx,t)-^{x,t)])) 
Jo Jn ^ 



du \ 

+ {l-C{xf)[u{x,t) -ij{x,t)]—]dxdt 



>0. 



The last integral is non-negative, because 

will do as a test function in the equation (7). This observation is important 
here. 

Aiming at difference quotients we give x the increment h. The translated 
function u{x + h,t) solves the obstacle problem with the translated obstacle 
iIj{x + h,t), all this with respect to the shifted domain Q^ x (0,T) where 
Q^ = {x\ X + h & Q}. For sufficiently small h we have 

/ / ({Ae{x + h,t),Vri{x,t))+r]{x,t)—^^^^^—^)dxdt>0 (8) 
Jo Jnh ^ ot J 

whenever ri{x,t) > ip{x + h,t) — u{x + h,t) and r] = on the parabolic 
boundary of Q^ x (0,T). Here 

Tj = ip{x + h,t) — u{x + h,t) + ({xY[u{x, t) — ip{x, t)] 



= C{xf [u{x + h,t)- u{x, t)] - (P{x) [ipix + h,t)- ■^{x, t)] 
-(1 - C{xY) [u{x + h,t)- ij{x + h, t)] 

will do. We obtain 

~l I ^ ((A.(a: + h,t), V(CAhw)) + CA^^u ^!^(^±M j dx dt 
+ 11^ ((A.(a: + h,t), V(CAh^)) + eA^^p du{x + h,t) ^ ^^ ^^ 

>[[ ({Mx + h,t),V{{l-axrMx + h,t)-ij{x + h,t)])) 
Jo Jnf^ ^ 

du{x + h,t)' 

> 0. 



+ (1 - C{x)P) [u{x + h,t)-^{x + h, t)] M^±M) j dx dt 



The last integral is positive because 

(1 - C{xY) Hx + h,t)- ij{x + h, t)] 

will do as a test function in the translated equation (8). This observation 
is essential here. The integrals in the left-hand member of the inequality 



are, in fact, taken only over the support of the function ({x). Hence we 
have an inequahty with integrals taken only over Qt, provided that \h\ < 
dist(suppC, i9f2). Thus Q^ is no longer directly involved. 

We add the two estimates, grouping the differences, and obtain 



+ / {Ae{x + h,t)-Ae{x,t),V{CAi,u))dxdt 
Jo Jn 

< / {Ae{x + h,t)- Ae{x,t),V{CAh'4j))dxdt 
Jo Jn 




CAY,U-AY,(^-^^dxdt 




C^Ah^-Ah(^)rfxrft. 



'0 Jn xL/t / Jo Jn 

The integrals with the time derivatives can be integrated by parts: 

^pd{A^ ^^ rft + /^ / C^AhV^ ■ Ah i^] dx dt 




Jn 



dt 2 

(Ahw)^ 



Jn 



'\dtJ 






dx + / C(x)Ah'ilj ■ AhM 
Jn 



dx 




eA^u-AJ^)dxdt. 
Jn ^dtJ 

Since AhW = Ah'?/' when t = 0, the above expression is majorized by 

i f C^((Ah^)^ - (Ah^)^) dx + ^ f I C{{^^.nf + (Ah^)') dxdt, 

where the inequality 2A\{a At^ip < (Ahw)^ + {Ay^ipY was used at time T. 

At this stage there are no "forbidden" time derivatives left and so we 
may safely let e go to zero. By Lemma 3 we may pass to the limit under the 
integral sign and hence the estimate for the limit u (no longer u^) becomes 




((AhA,V(C^AhM))rfxrft 
'0 Jn 

rp 

< [ /((AhA,V(C^Ah^))rfxrft 
Jo Jn 



+ 



\j j C (( Ah«)' + ( Ah^) ')dxdt+]^j C^( Ah^)| dx, 



where 

AhA = A(a; + h,t) - A{x, t). 

We write this more conveniently as 

C^((AhA,VAhM))rfxdt 




Jn 



I 



< 



+ 



/ [ pC'^\AhA\\Ahu\\VC\dxdt 
Jo Jn 



II 



[ [ pC'^\Ai,A\\Ai,^\\VC\dxdt 
Jo Jn 



(9) 



III 



+ f f C\^hA\\VAhip\dxdt 
Jo Jn 

+ : 



\j j C^((Ah«)2(Ah^)') dxdt+\j Ci^^^frdx. 
The integrand on left-hand side is (AhA, VAhw)) = 



, P-2 . 



, p-2 . 



{\Wu{x + h,t)\-^Wu{x + h,t) - \Wu{x,t)\-^Wu{x,t),Wu{x+h,t) -Vu{x,t)) 

(10) 



> 4 |F(x + /i, t) - F(x, t)|' = 4l^hFr, 



where the elementary inequality 



pz 



, p-2. 



, p-2 



|6| Vfe - |a| Va < {\b\P-% - \a\P-^a, b - a) 



for vectors was used. We aim at an estimate for the integral of C^| AhF|. 

We divide the Ah-terms by \h\ so that the desired difference quotients 
appear. The estimate 



AuA 



h 



< (P - 1) 



AhF 



h 



\Vu{x + h,t)\^ + \Vu{x,t)\^), 



coming from the elementary vector inequality 



\\bf-%-\a\P-^a\ < {p-l)(\b\^ + \a\'^) {bl^b-lal'^a 
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is used in the integrands of I, II, and III. In I we split the factors so that 



< p{p - 1) 



h ' 



h ' 



A.A 



h 



AhU 



IVCI 



and use Young's inequality 



IVCI 



C"^(|VM(x,t)|^ 



2p 



p^^ 



\Vu{x+h,t)\—) 



abc < 



e^a^ e-PbP {p - 2)c~^ 



p 



to get the bound 



\h\ 



I ^p{p-iy f^ I ^, 




Jn 



2p 



AhF 



h 



dxdt 



+ {P- 1)^'*' 




Jn 



AhM 



\VC\^dxdt 



+Cb 



I f C' (l Vm(x, t) f + \Vu{x + h, t) \A dx dt 



The integral II/|/ip has a similar majorant, the only difference being that 
AhM be replaced by Ah'?/'. The integrand of III is estimated in a similar way: 



< p{p - 1) 



■^£|AhF 
(2I 



< 



' h 
p{p - l)e 



CV 



■e 



pC 

AhM 

^ h 
AhF 



A.A 



h 



Ah^ 



h 



_p-2 



, p-2 



, p-2^ 



+ {p- l)e-X- 



+CpC{\'^uix,t)f + \Vuix + h,t)\P). 



C~{\Vu{x,t)\^- + \Vu{x + h,t)\ — ) 
AhV" " 



h 



11 



Adding up the three integrated estimates, we arrive at 
I + 11 + III 




Jn 

T 



< 



AhF 



dx dt 



+ '^-i'^"/„/„( 



Jn 



AhM 



h 



ivcr+ 



AhV^ 



h 



ivcr+c^ 



V(AhV^ 



dxdt 



+ 



3cp f f c(\Vu{x,t)\P+\Vu{x + h,t)\p)dxdt. 



This complements (9). Recall (10). The next step is to absorb the first 
integral above in the right-hand member into the minorant in (10) by fixing 
e small enough, say 

^ p{p-l)e' _ 2 
2 p2- 

The resulting estimate, written out without abbreviations, is 



^ ' F{x + h,t)-F{x,t) 




Jn 



<ar. 



h 



dxdt 



T 




'0 Jn 



u{x + h,t) — u{x,t) 



+a. 




h 
ip{x + h,t) — ip{x,t) 



'0 Jn 

+ap lie 
Jo Jn 

+' 



h 



\V CI" dxdt 



\VC\^dxdt 



Vip{x + h,t) -V^{x,t) 



h 



dxdt 



K [ I C(\^u{x,t)f+\Vu{x + h,t)\Adxdt 




+cp I I e 



'0 Jn 

i-T 



u{x + h,t) — u{x,t) 



+ Cri 




e 



Jn 



h 
iptix + h,t) -ilJt{x,t) 



h 



dxdt 

2 

dxdt 



+Cp / C^ 
'n 



ipi^x + h,T) — ip{x, T) 



h 



dx. 



12 



where the constants depend only on p. Finally, letting the increment h ^ 
in any desired direction, we arrive at the estimate in the theorem. Here we 
use the characterization of Sobolev spaces in terms of integrated differential 
quotients, cf. [G, Chapter 8.1]. This concludes our proof of Theorem 5. 

Corollary 6 Ifuis the solution to the obstacle problem with the obstacle ip, 
then ApU = V ■ (|V'u|^~^V'u) belongs to the space L^^^ (^t) CLnd 

nip ApU dxdt= / ( \Vu\P~'^Vu, Vv?) dx dt 
! Jo Jn 

for all test functions ip in C^{Qt)- 

Proof: Since F is in Sobolev's space and p /, 2, we can differentiate 



£-2. 



Vu\P-^Vu= |F|~F 



and hence 



d 



dx. 



{\Vu\P-^Vu) 



< 2(1- - F ~ 
p/ 



dF 



dx^i 



By Holder's inequality 



dXn 



(|v«r^v«)GL-;(fiT) 



since F G /.^(f]^) and DY e L'^inr). □ 



4 The Time Derivative 

For the proof of the Theorem we notice that the contact set S = {u = ip} is 
a closed subset of Qt and that its complement T = ^2^ \ S is open. In the 
set T, where the obstacle does not hinder, m is a solution to the Evolutionary 
p-Laplace Equation ut = ApU. In other words, whenever (p G C^(T), 

nu—dxdt= I I ( |Vm|^^^Vm, V0) (ixiit = — / / (pApudxdt, 
! ^^ Jo Jn Jo Jn 

the actual set of integration being T. Here Corollary was used. Thus Ut is 
available, but only in T to begin with. (See also [L2].). Let denote an 
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arbitrary test function in C^{Qt)- We need a specific test function with 
compact support in T. To construct it, define 

9k = min{l, k{u — ip)}, /c = 1, 2, ■ ■ ■ . 

Then 1 — ^^ = 1 in S and pointwise the monotone convergence 1 — Ok — > X'b 
holds. Moreover, the support of 6k is compact in T. The time derivative of 
6k is available! 
Using 

[ I u—{6k(t))dxdt= I [ (\Vu\P-^Vu,V(6k(l)))dxdt, 
Jo Jn dt Jo J^ 

we write 

n0ApM dxdt = - { I Vm|^" Vm, V0) dx dt 

_! Jo Jn 



- / f {\VuY'-^Vu,V{6k<p+{l-6k)<p))dxdt 
Jo Jn 



f f {\Vu\P-^Vu,Vi6k(j)))dxdt- f f {\Vu\P-^Vu,V{{l-6k)(l)))dxdt 
'0 Jn Jo Jn 



nu—{6k4>) dxdt+ / (1 - 6k)(j)ApU dx dt. 
^ A ^^ Jo Jn 



The last integral has the limit 



T 




lim / / (I — 6k)(p Apudx dt = / (p Ap^p dx dt. 



^=^0 -'0 Jn 



In the integral with the time derivative we write 

-u-{6k<P) = -u^ + {u- ^)-(l - 6k)<P + ^^(1 - 6k) 



and obtain 



T r a pT 




u-—{6k4>)dxdt= / / —u-^dxdt 
'0 Jn dt Jo J^ dt 

+ f f{u-^)^iil-6k)(l))dxdt-f fil-6k)(j)^dxdt, 
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where an integration by parts has produced the last integral It has the 
evident hmit 

f'^ f dip f f dip 

hm / / {\ — 6k)4>-;T- dxdt = I I (p-— dxdt. 



k^o Jq J^ dt J J dt 



The middle integral vanishes as A; — )■ 0: 

r f, ^.d 



Jo Jn 


^'dt^' 


K-^ "K/ 


'T/ ""^ "" 












Jo Jn 


-ip){i- 


- ^'^ dt 


dxdt — 


f 

Jo . 


/ (piu — ih)—:— dxdt 
In dt 




77<- 

Jo Jn 


-ip){i- 


" ^'^ dt 


dxdt — 


1 
2k ^ 


11 

'0 Jn 


<p—dldxdt 






77<" 

Jo Jn 


-ip){i- 


- ^'^ dt 


dxdt + 


1 
2k ^ 


11 

'0 Jn 


el ^j" dxdt 

'^ dt 


— y 

k—¥oo 


+ 


Collecting 


; results, 
















l-T f. 






fT f. 






r r 






/ / 


Ap-u dxdt = - 


/ / 


(pt dx dt — 


i^t-A, 


■jip)(pdxdt. 


Jo Jn 






Jo Jn 






J J 







In other words, the final formula 

— / / u(pt dxdt = I (p [Ap-u + {ipt — Apip)x~\ dx dt 
Jo Jn Jo Jn 

holds for every (p in C^{VLt)- Therefore 

ut = ApU + {ipt - Ap^)xs 

and this is a function belonging to L^^^^ (^t)- This concludes the proof of 
Theorem 2. D 

Peter Lindqvist 
Department of Mathematical Sciences, Norwegian University of Science and Techno- 
logy, N-7491, Trondheim, Norway 
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